Heavy metal concentration according to shrimp species and organ specificity: Monitoring and human risk assessment.

This study assessed heavy metal levels (lead (Pb), cadmium (Cd), total arsenic (tAs), arsenite (As (III)), arsenate (As (V)), monomethyl arsenic acid (MMA), dimethylarsinic acid (DMA), total mercury (tHg), and methylmercury (MeHg)) in six organs (total portion, head, body, shell, muscle, and intestine) of 11 shrimp species distributed in Korea. Shrimp exhibited significant variability in heavy metal accumulation, with Alaskan pink and dried shrimp (Lesser glass, Southern rough, and Chinese ditch prawn) showing the highest metal concentrations. Notably, the intestine having the highest overall metal content, while Cd was most prominent in the head, tHg was highest in the muscle. The Hazard Quotient values of 11 shrimp species in South Korea were below the European Food Safety Authority's allowable limits for heavy metals. This study illuminates the heavy metal profiles of distributed shrimp in Korea and emphasizes the ongoing need for monitoring heavy metals on seafood to ensure consumer safety.

Monitoring arsenic species concentration in rice-based processed products distributed in South Korean markets and related risk assessment.

Rice is an important grain as a major source of carbohydrates in Asia but contains more arsenic (As) than other grains. A total of 239 rice-based processed foods (rice, n = 30; rice cake, n = 30; porridge, n = 39; noodles, n = 33; bread, n = 20; snack, n = 59; powder, n = 28) were purchased in 2019 from domestic markets to measure total As (tAs) and As species. The average tAs and inorganic As (iAs) in each sample group ranged from 20 to 180 µg/kg (porridge for baby to noodle) and 4.4-85 µg/kg (porridge for baby to powder), respectively. The correlation between the iAs and tAs was affected by the variety of ingredients, such as the presence of seaweed (tAs) and the milling type of rice (iAs). Although rice cakes and baby rice-based powders are a source of concern for both adults and children, respectively, risk assessments indicate that most rice-based foods are generally safe to consume in South Korea. Supplementary Information: The online version contains supplementary material available at 10.1007/s10068-023-01270-9.

Risk assessment of polycyclic aromatic hydrocarbons in meat and edible oils: results of a total diet study in South Korea.

Polycyclic aromatic hydrocarbons (PAHs) constitute carcinogens. In this study, the risk of PAHs being consumed through meat and edible oils was assessed using a total diet study. Results were monitored by applying the toxic equivalency factor of benzo[a]pyrene; among each category, this factor was highest in grilled beef chitterlings (1.35 µg/kg), grilled Wiener sausages (1.20 µg/kg), fried chicken wings (0.70 µg/kg), and stir-fried perilla oil (1.29 µg/kg). The chronic daily intake was calculated, and risk characterization was estimated by applying the margin of exposure using the benchmark dose approach. Most samples analyzed in our study were denoted as having no concern; however, the intake group of stir-fried beef chitterlings, pan-fried pink sausage, deep-fried pork loin, and grilled duck was regarded as possible concern, and grilled chicken was assessed as having low concern. PAH changes must be monitored on a regular basis. Supplementary Information: The online version contains supplementary material available at 10.1007/s10068-022-01137-5.

Method Validation for Determination of Thallium by Inductively Coupled Plasma Mass Spectrometry and Monitoring of Various Foods in South Korea.

Thallium (Tl) is a rare element and one of the most harmful metals. This study validated an analytical method for determining Tl in foods by inductively coupled plasma mass spectrometry (ICP-MS) based on food matrices and calories. For six representative foods, the method's correlation coefficient (R2) was above 0.999, and the method limit of detection (MLOD) was 0.0070-0.0498 µg kg-1, with accuracy ranging from 82.06% to 119.81% and precision within 10%. We investigated 304 various foods in the South Korean market, including agricultural, fishery, livestock, and processed foods. Tl above the MLOD level was detected in 148 samples and was less than 10 µg kg-1 in 98% of the samples. Comparing the Tl concentrations among food groups revealed that fisheries and animal products had higher Tl contents than cereals and vegetables. Tl exposure via food intake did not exceed the health guidance level.

Polycyclic aromatic hydrocarbons in beverage and dairy products in South Korea: a risk characterization using the total diet study.

Polycyclic aromatic hydrocarbons (PAHs) were analyzed using gas chromatography-mass spectrometry in 115 dairy products and beverages, including alcoholic, grain, carbonated, and functional drinks; fruit and vegetable juices; coffee; and tea, purchased from 10 local city markets in South Korea. The sample groups were divided into non-fatty and fatty groups, pretreated with the ultrasound-assisted extraction method and saponification method, respectively. The limit of detection, limit of quantification, and accuracy were 0.038-0.185 µg/kg, 0.114-0.560 µg/kg, and 87.64-112.25%, respectively. The measurement uncertainty was ≤ 6.38% for eight PAHs (PAH8). PAH8 was detected in 41 of the 115 samples, ranging from 0.041 to 7.793 µg/kg. The risk assessment revealed that the margin of exposure for PAH8 ranged from 3.60 × 104 to 7.84 × 1011 in the mean intake groups and from 3.60 × 104 to 5.33 × 1011 in the P97.5 intake groups. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10068-021-00927-7.

Origin and organ-specific bioaccumulation pattern of perfluorinated alkyl substances in crabs.

Consumption of seafood is a major contributor to perfluorinated alkyl substances (PFASs) exposure. Crabs contain high levels of PFASs, and different PFASs are concentrated in their tissues depending on their habitat. Despite South Korea importing huge quantities of crabs, no investigation has been conducted on the effect of PFAS exposure. This study investigated the risk of exposure to PFASs when ingesting crabs. To determine the risk of exposure, 19 different PFAS species were measured in the edible parts (body, legs, offal, and eggs) of crabs originating from South Korea (n = 17), China (n = 14), India (n = 7), and Pakistan (n = 31), which were distributed in the fish markets of South Korea. The results revealed that, in contrast to short-chain PFASs, long-chain PFASs (PFCAs≥8, PFSAs≥6, and perfluorooactane sulfonamidoacetic acids (FOSAAs)≥8) were detected in crab samples from all four countries of origin, and in all the edible parts except for the legs. Perfluorooctanoic acid (PFOA; 16.9 ng/g in South Korea, 9.42 ng/g in China) and perfluoro-n-tridecanoic acid (PFTrDA; 5.35 ng/g in South Korea, 2.40 ng/g in China) were the predominant perfluoroalkyl carboxylic acids (PFCAs) detected in the crabs originating from South Korea and China, and perfluorooctane sulfonic acid (PFOS; 7.02 ng/g in Pakistan, 5.88 ng/g in India) was the predominant perfluoroalkyl sulfonic acids (PFSAs) detected in crabs originating from Pakistan and India. These results indicate that PFASs that are accumulated in crabs differ depending on the ocean from which they originate. The concentrations of PFOA and PFOS were significantly higher in the eggs and offal than in the legs and body of the crab. The average daily intake of PFOA and PFOS in South Koreans ranges from 0.01% to 0.07% based on the tolerable daily intake of EFSA and MFDS. These results establish the PFAS profiles and risk assessment of crabs that are distributed in South Korea.

Comparing levels of perfluorinated compounds in processed marine products.

Perfluorinated compounds (PFCs) are widely distributed in nature and have many applications due to their unique chemical and physicochemical properties. While, PFCs are present in soil, water, and air, their pathway for entry into the human body is circumstantially via contaminated food. The demand for seafood has been increasing. In this study, we investigated via LC-MS/MS, the content of 19 different types of PFCs in 302 samples belonging to five different categories of the typical South Korean seafood. The highest levels of PFOA, PFTrDA, PFOS, and PFPeA were found in dried seafood, canned and salted seafood, processed fish food, and seasoned laver, respectively. The levels of PFOA and PFOS were compared on the basis of various criteria including the nomenclature, biological classification, and habitat of the source of seafood. High levels of PFOA and PFOS were detected in anchovy, sea squirt, and mackerel based on the nomenclature of raw seafood, in crustaceans based on the biological classification, and in demersal organisms based on the habitat. The human intake values of PFOA and PFOS from the processed marine products in South Korea were lower than the tolerable daily intake, indicating that the consumption of these processed seafood poses no immediate harm.

Identification of environmental determinants for spatio-temporal patterns of norovirus outbreaks in Korea using a geographic information system and binary response models.

Although norovirus outbreaks are well-recognized to have strong winter seasonality relevant to low temperature and humidity, the role of artificial human-made features within geographical areas in norovirus outbreaks has rarely been studied. The aim of this study is to assess the natural and human-made environmental factors favoring the occurrence of norovirus outbreaks using nationwide surveillance data. We used a geographic information system and binary response models to examine whether the norovirus outbreaks are spatially patterned and whether these patterns are associated with specific environmental variables including service levels of water supply and sanitation systems and land-use types. The results showed that small-scale low-tech local sewage treatment plants and winter sports areas were statistically significant factors favoring norovirus outbreaks. Compactness of the land development also affected the occurrence of norovirus outbreaks; transportation, water, and forest land-uses were less favored for effective transmission of norovirus, while commercial areas were associated with an increased rate of norovirus outbreaks. We observed associations of norovirus outbreaks with various outcomes of human activities, including discharge of poorly treated sewage, overcrowding of people during winter season, and compactness of land development, which might help prioritize target regions and strategies for the management of norovirus outbreaks.

Deep hypothermic circulatory arrest versus non-deep hypothermic beating heart strategy in descending thoracic or thoracoabdominal aortic surgery.

OBJECTIVES: The ideal cardiopulmonary bypass (CPB) strategy during open surgical repair of the descending thoracic aorta (DTA) or thoracoabdominal aorta (TAA) is controversial. This study aimed to compare the clinical outcomes between deep hypothermic circulatory arrest (DHCA) and non-deep hypothermic beating heart CPB (non-DHCA) for DTA or TAA replacement. METHODS: From January 1994 to August 2011, 259 patients underwent DTA or TAA replacement. Of these, 212, who were judged to be suitable for both DHCA (n = 79) and non-DHCA (n = 109), were analysed. In-hospital outcomes were compared using propensity scores and inverse-probability-weighting adjustment based on 20 preoperative variables to reduce treatment selection bias. RESULTS: Early mortality was 12.7% in the DHCA group and 7.5% in the non-DHCA group (P = 0.23). Major adverse outcomes included stroke in 13 patients (6.1%), paraplegia in 10 (4.7%), low cardiac output syndrome (LCOS) in 17 (8.0%) and multiorgan failure in 12 (5.7%). After adjustment, patients who underwent DHCA were at a risk of death (odds ratio (OR), 1.86; P = 0.18) and permanent neurological injury (OR, 1.06; P = 0.90) similar to that of those who underwent non-DHCA, but at greater risk of LCOS (OR, 3.85; P = 0.012). Furthermore, prolonged ventilator support (>24 h) was more frequent with DHCA than with non-DHCA (OR, 2.33; P = 0.004). CONCLUSIONS: Compared with non-DHCA, DHCA was associated with greater risk of postoperative LCOS and prolonged ventilator support. Therefore, non-DHCA seems to be a more appropriate option than DHCA for open DTA/TAA repair whenever the aortic anatomy lends itself to this approach.

Flexible marginalized models for bivariate longitudinal ordinal data.

Random effects models are commonly used to analyze longitudinal categorical data. Marginalized random effects models are a class of models that permit direct estimation of marginal mean parameters and characterize serial correlation for longitudinal categorical data via random effects (Heagerty, 1999). Marginally specified logistic-normal models for longitudinal binary data. Biometrics 55, 688-698; Lee and Daniels, 2008. Marginalized models for longitudinal ordinal data with application to quality of life studies. Statistics in Medicine 27, 4359-4380). In this paper, we propose a Kronecker product (KP) covariance structure to capture the correlation between processes at a given time and the correlation within a process over time (serial correlation) for bivariate longitudinal ordinal data. For the latter, we consider a more general class of models than standard (first-order) autoregressive correlation models, by re-parameterizing the correlation matrix using partial autocorrelations (Daniels and Pourahmadi, 2009). Modeling covariance matrices via partial autocorrelations. Journal of Multivariate Analysis 100, 2352-2363). We assess the reasonableness of the KP structure with a score test. A maximum marginal likelihood estimation method is proposed utilizing a quasi-Newton algorithm with quasi-Monte Carlo integration of the random effects. We examine the effects of demographic factors on metabolic syndrome and C-reactive protein using the proposed models.

Traffic violations versus driving errors of older adults: informing clinical practice.

Certain driving errors are predictive of crashes, but whether the type of errors evaluated during on-road assessment is similar to traffic violations that are associated with crashes is unknown. Using the crash data of 5,345 older drivers and expert reviewers, we constructed a violation-to-error classification based on rater agreement. We examined the effects of predictor variables on crash-related injuries by risk probability using logistic regression. Drivers' mean age was 76.08 (standard deviation = 7.10); 45.7% were women. Of drivers, 44.6% sustained crash-related injuries, and female drivers had a higher injury probability (44%) than male drivers (29%). Lane maintenance, yielding, and gap acceptance errors predicted crash-related injuries with almost 50% probability; speed regulation (34%), vehicle positioning (25%), and adjustment-to-stimuli (21%) errors predicted crash-related injuries to a lesser degree. We suggest injury prevention strategies for clinicians and researchers to consider for older drivers, especially older women.

Validity and critical driving errors of on-road assessment for older drivers.

OBJECTIVES: We examined the validity of our on-road driving assessment to quantify its outcomes. METHOD: Older drivers (N = 127) completed a driving assessment on a standardized road course. Measurements included demographics, driving errors, and driving test outcomes; a categorical global rating score (pass-fail); and the sum of maneuvers (SMS) score (0-273). RESULTS: There were significant differences in the SMS (F = 29.9, df = 1, p < or = .001) between drivers who passed the driving test and those who failed. The SMS cutoff value of 230 points was established as the criterion because it yielded the most optimal combination of sensitivity (0.91) and specificity (0.87). The strongest predictors of failure were adjustment to stimuli and lane maintenance errors. CONCLUSION: The SMS differentiated between passing and failing drivers and can be used to inform clinical decision making.

Bayesian model-based tight clustering for time course data.

Cluster analysis has been widely used to explore thousands of gene expressions from microarray analysis and identify a small number of similar genes (objects) for further detailed biological investigation. However, most clustering algorithms tend to identify loose clusters with too many genes. In this paper, we propose a Bayesian tight clustering method for time course gene expression data, which selects a small number of closely-related genes and constructs tight clusters only with these closely-related genes.

Hydrochemical and multivariate statistical interpretations of spatial controls of nitrate concentrations in a shallow alluvial aquifer around oxbow lakes (Osong area, central Korea).

Hydrochemical and multivariate statistical interpretations of 16 physicochemical parameters of 45 groundwater samples from a riverside alluvial aquifer underneath an agricultural area in Osong, central Korea, were performed in this study to understand the spatial controls of nitrate concentrations in terms of biogeochemical processes occurring near oxbow lakes within a fluvial plain. Nitrate concentrations in groundwater showed a large variability from 0.1 to 190.6 mg/L (mean=35.0 mg/L) with significantly lower values near oxbow lakes. The evaluation of hydrochemical data indicated that the groundwater chemistry (especially, degree of nitrate contamination) is mainly controlled by two competing processes: 1) agricultural contamination and 2) redox processes. In addition, results of factorial kriging, consisting of two steps (i.e., co-regionalization and factor analysis), reliably showed a spatial control of the concentrations of nitrate and other redox-sensitive species; in particular, significant denitrification was observed restrictedly near oxbow lakes. The results of this study indicate that sub-oxic conditions in an alluvial groundwater system are developed geologically and geochemically in and near oxbow lakes, which can effectively enhance the natural attenuation of nitrate before the groundwater discharges to nearby streams. This study also demonstrates the usefulness of multivariate statistical analysis in groundwater study as a supplementary tool for interpretation of complex hydrochemical data sets.

Marginalized random effects models for multivariate longitudinal binary data.

Generalized linear models with random effects are often used to explain the serial dependence of longitudinal categorical data. Marginalized random effects models (MREMs) permit likelihood-based estimations of marginal mean parameters and also explain the serial dependence of longitudinal data. In this paper, we extend the MREM to accommodate multivariate longitudinal binary data using a new covariance matrix with a Kronecker decomposition, which easily explains both the serial dependence and time-specific response correlation. A maximum marginal likelihood estimation is proposed utilizing a quasi-Newton algorithm with quasi-Monte Carlo integration of the random effects. Our approach is applied to analyze metabolic syndrome data from the Korean Genomic Epidemiology Study for Korean adults.

Clinical predictors of older driver performance on a standardized road test.

OBJECTIVES: To determine the relationship between clinical variables (demographics, cognitive testing, comorbidities, and medications) and failing a standardized road test in older adults. METHODS: Analysis of on-the road studies performed in optimal weather conditions, between January 1, 2005, and May 1, 2007. The standardized testing was held at the National Older Driver Research and Training Center (NODRTC), Florida, and included 127 adults aged 65 and older with current driver licenses, recruited by advertisement from the Gainesville, Florida, community. Measurements consist of demographics, self-reported medications and medical conditions, cognitive testing including Trail Making Part B, global rating score (pass/fail), and driver maneuver score (0-273, with 273 indicating perfect driving or zero errors). RESULTS: A total of 127 older adults completed the protocol. Mean age was 74.8 years (SD = 6.3); 46.5% females. Mean time for Trail Making Part B was 114.3 seconds (SD of 83). Among the 127 drivers, the mean Sum of Maneuvers Score was 238.9 (SD of 25.0) and 24 (19%) failed the driver test. Odds ratio estimates for failing the test included advanced age (6.7, 95% CI 2.2 to 19.8), presence of a neurological disease (2.8, 95% CI 1.2 to 6.5), and prolonged time to complete the Trail Making Part B cognitive test (2.5, 95% CI 1.0 to 5.9). Conversely, odds ratio estimates lowering the risk of failure included taking a non-diabetic hormonal medications (e.g., thyroid and estrogen drugs; 0.3, 95% CI .09 to 0.7) and having a musculoskeletal diagnosis (0.3, 95% CI .1 to 0.7). CONCLUSIONS: To our knowledge, this is the first study to examine the medical predictors of failing a standardized road test. Advanced age and prolonged time on Trail Making Part B were the two major predictors of test failure and a lower Sum of Maneuvers Score. Our study also found that having a neurological diagnosis (primarily cerebrovascular and Parkinson's disease) predicted test failure. Medications from neurological class also predicted a lower Sum of Maneuvers Score. Further study needs to be done to explain the apparent protective effect of musculoskeletal conditions and hormonal medications.

Model-based Bayesian clustering (MBBC).

MOTIVATION: The program MBBC 2.0 clusters time-course microarray data using a Bayesian product partition model. RESULTS: The Bayesian product partition model in Booth et al. (2007) simultaneously searches for the optimal number of clusters, and assigns cluster memberships based on temporal changes of gene expressions. MBBC 2.0 to makes this method easily available for statisticians and scientists, and is built with three free computer language software packages: Ox, R and C++, taking advantage of the strengths of each language. Within MBBC, the search algorithm is implemented with Ox and resulting graphs are drawn with R. A user-friendly graphical interface is built with C++ to run the Ox and R programs internally. Thus, MBBC users are not required to know how to use Ox, R or C++, but they must be pre-installed. AVAILABILITY: A self-extractable zip file, MBBC20zip.exe, is available at the MBBC webpage www.stat.ufl.edu/~casella/mbbc/, which contains MBBC.exe, source files, and all other related files. The current version works only in the Windows operating system. A free installation program and overview for Ox is available at www.doornik.com. A detailed installation guide for Ox is provided by MBBC, and is accessible without installing Ox. R is available at www.r-project.org/.

Phosphorus and other soil components in a dairy effluent sprayfield within the central Florida Ridge.

There is concern that P from dairy effluent sprayfields will leach into groundwater beneath Suwannee River basins in northern Florida. Our purpose was to describe the effects of dairy effluent irrigation on the movement of soil P and other nutrients within the upper soil profile of a sprayfield over three 12-mo cycles (April 1998-March 2001). Effluent P rates of 70, 110, and 165 kg ha(-1) cycle(-1) were applied to forages that were grown year-round. The soil is a deep, excessively drained sand (thermic, uncoated Typic Quartzipsamment). Mean P concentration in soil water below the rooting zone (152-cm depth) was < or = 0.1 mg L(-1) during 11 3-mo periods. Mehlich-1-extractable (M1) P, Al, and Ca in the topsoil increased over time but did not change in subsoil depths of 25 to 51, 51 to 71, 71 to 97, and 97 to 122 cm. Topsoil Ca increased as effluent rate increased. High Ca levels were found in dairy effluent (avg.: 305 mg L(-1)) and supplemental irrigation water (avg.: 145 mg L(-1)) which likely played a role in retaining P in the topsoil. An effect of effluent rate on P and Al concentrations in the topsoil was not detected, probably due to large and variable quantities present at project initiation. The P retention capacity (i.e., Al plus Fe) increased in the topsoil because Al increased. Dairy effluent contained Al (avg.: 31 mg L(-1)). Phosphorus saturation ratio (PSR) increased over time in the topsoil but not in subsoil layers. Regardless of effluent rate, the P retention capacity and PSR of subsoil, which contained 119 to 229 mg kg(-1) of Al, should be taken into account when assessing the risk of P moving below the rooting zone of most forage crops.

Normalization of dye bias in microarray data using the mixture of splines model.

This paper discusses characteristics of dye biases in microarray data that the conventional normalization methods do not handle, and proposes a new normalization method involving a mixture of splines model. We also develop a test for between-group comparisons of each gene that is designed to be used with our proposed method.

Five year-round forage systems in a dairy effluent sprayfield: phosphorus removal.

In northern Florida, forages are grown in dairy effluent sprayfields to recover excess P. Our purpose was to evaluate five year-round forage systems for their capacity to remove P from a dairy sprayfield. The soil is a Kershaw sand (thermic, uncoated Typic Quartzipsamment). Systems included bermudagrass (Cynodon spp.)-rye (Secale cereale L.) (BR), perennial peanut (Arachis glabrata Benth.)-rye (PR), corn (Zea mays L.)-forage sorghum [Sorghum bicolor (L.) Moench]-rye (CSR), corn-bermudagrass-rye (CBR), and corn-perennial peanut-rye (CPR). Forages were grown for five 12-mo cycles. Effluent P rates were 80, 120, and 165 kg ha-1 cycle-1. The 5-cycle P removal was 67 kg ha-1 cycle-1 for BR, 54 kg ha-1 for CBR, 52 kg for CSR, 45 kg for PR, and 43 for CPR. Removal of P by winter rye was low. There were differences in system rankings among cycles primarily due to changes in the performance of perennial forages. In the first two cycles, BR had the greatest P removal (91 kg ha-1 cycle-1) due to high bermudagrass yield and P concentration. In the first cycle, P removal was lowest for PR (36 kg ha-1) because perennial peanut was slow to establish. In later cycles, P removal for BR declined because bermudagrass yield and P concentration declined. It increased for PR because peanut yield increased. The yield of corn in CBR, CPR, and CSR was consistently high but P concentration was modest (avg. 2.2 g kg-1). Sorghum produced moderate but stable yield and had low P levels (avg. 1.8 g kg-1). Effluent rate marginally affected the performance of most grasses. For P recovery in dairy sprayfields in northern Florida, the best warm-season forage would likely be a high yielding, persistent bermudagrass.